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Pick a basis to express an arbitrary tensor, then form a polynomial in its coefficients. 
Any action of a group on the tensor induces an action of this group on the polynomial. 
A polynomial unchanged under the group's action is termed a polynomial invariant. The 
theory is rich and the study of invariants dates back several centuries. The topic has recently 
attracted attention also in the study of quantum theory 

Finding a complete and minimal set of polynomial invariants of the local unitary group 
for a quantum state would allow one to write it in a reduced coordinate system, without ref- 
erence to local degrees of freedom. This basis independent description aids in understanding 
properties that are invariant under the group action, including entanglement measures and 
entropies. While these invariants could deliver new structural insights into the underlying 
properties of non-locality of a quantum state, the polynomial equations that describe them 
contain a large number of terms with coefficients raised to high powers. The complexity of 
the description of the invariants themselves limits their usefulness, their physical interpreta- 
tions and mathematical progress in the field. We seek to address this problem by expressing 
the invariants in terms of tensor networks. 

In this paper we present a variant of the graphical tensor calculus of Penrose 0] for 
the purpose of representing and computing polynomial invariants of quantum states and 
tensor networks in general. The underlying mathematical structure of the physics being 
described by the invariants is reflected in the structure of the resulting tensor networks. 
By using specific graphical rewrite rules (tensor contraction rules), our methods enable one 
to contract and simplify a tensor network representing any given polynomial invariant of 
a bipartite pure state to the point where the network is succinctly expressed in terms of 
Schmidt coefficients. This serves as a graphical proof of the invariance of the quantity 
represented by the network, as well as a conceptual aid geared towards understanding the 
meaning behind the invariants. 

The proposed method seems to open several avenues related to both invariant theory and 
tensor network states. We intend to make the present paper serve as the backbone of several 
results related to these developments. It contains the framework, key definitions and several 
results related to matrix product states (MPS) 0, IsJ. At the heart of MPS is the tensor 
network description of repeated bipartitions of a quantum state. By capturing the singular 
value decomposition in a tensor network where all internal components have clearly defined 
algebraic properties, we present some improvements in the graphical tensor calculus used to 
describe matrix product states, as well as invariants in general. In this regard, our results 
on matrix product states take an important first step in uniting invariant theory with tensor 
network states in a more general setting. 

An important future goal is the development of a theory of efficiently calculable invariants 
for tensor network states in geometries more complicated than ID matrix product states. 
The methods we present here are proven to work for matrix product factorizations of quan- 
tum states into bipartitions, but the problem becomes more difficult in the general setting. 
In this sense, here we are communicating partial results towards the larger goal of developing 
a full and unified invariant theory for tensor network states in general. 

One might use tensor network algorithms 0-0] to design and contract invariants of in- 
terest to condensed matter physics, suggesting some practical utility of our results. Our 
methods offer a better conceptual aid to understand how the numerical value of the invari- 
ant relates to properties present in the state in question. This could lend in the development 
and understanding of new quantities of interest to condensed matter physics as well as a 
direct connection to tensor network state computer algorithms, developed to evaluate tensor 
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contractions. 

The structure of the present paper is as follows. We begin by recalling the fundamental 
notions of the tensor calculus in Sections [TJ [2] and [31 This leads to the diagrammatic SVD, 
which is used in Section H] to factor a given quantum state into a matrix product state. We 
then connect the tensor calculus to polynomial invariants in Section [5j Before concluding, we 
also consider the application of invariants to calculate entropies and entanglement measures. 



1. PENROSE GRAPHICAL NOTATION FOR TENSOR NETWORKS 



Penrose graphical notation [6| is a diagrammatic notation for tensor networks. This 
notation is becoming well known inside the tensor network algorithms community. It can 
make the manipulation of complicated tensor networks much easier and more intuitive. 
Contributions on the topic we found influential can be found in In our previous work, 

we have adapted the graphical notation and surrounding methods to describe generalized 
quantum circuits [l5|, tensor network states 15 -17], open quantum systems 18, 19 1 as well 
as decidability in algorithms based on tensor contractions [20|. 

In the string diagram notation, a tensor is a graphical shape with a number of input 
legs (or "arms") pointing up, and output legs pointing down. Individual arms as well as 
individual legs each independently correspond to an index. For example, 



(a) 



* 1 



(b) 



T 



J 1 1 k 



diagram (a) represents the tensor ipi and (b) the tensor T l j k . A tensor with n indices up and 
m down is called a valence- (n, m) tensor and sometimes a valence-/? tensor for k = n + m. 
Often, to conserve space, the diagrams are rotated 90 degrees counterclockwise. In quantum 
physics parlance one introduces a computational basis and expands the tensors in it; in 
which case T l - k is understood not as abstract index notation but as the actual components 
of the tensor: 

ijk 

In practice there is little room for confusion however. 

There are three special "wire tensors" that play the role of the metric tensorQ They are 
given diagrammatically as 



(a) 



si 



(b) 



u 



(c) 



n 



The identity tensor (a) is used for index contraction by connecting the corresponding legs, 
and the cup (b) and cap (c) are metric tensors used for raising and lowering indices. Ex- 



We will always work in a flat Euclidean space, which renders the metric tensors trivial. 
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(a) 



i 




J 


A 




T 


k 








B 



rn 




FIG. 1. Illustration of the graphical notation, (a) Contraction of tensor T with tensors A and B 
amounts to joining indices: T' l - k A 3 l B k m . (b) Permutation of indices by crossing wires: A deg B aC f. 



panding them in the computational basis we obtain 

ij k 

(ui = E^^'i = E^i' and 

ij k 



(2) 
(3) 
(4) 



ij 



(i) One can raise and subsequently lower an index or vice versa, which amounts essentially 
to doing nothing at all. This scenario is captured diagrammatically by the so called 

snake or zig-zag equation 



C 



together with its mirror image [6[. In tensor index notation, it is expressed succinctly 
as 8 ij 8 jk = 5\. 

(ii) Crossing two wires (as in diagram (a) below) is equivalent to swapping the relative 
order of the corresponding vector spaces. 



(a) 



X 



(b) 



XX = 



(b) illustrates that the swap operation is self inverse. It may be written as SWAP 1,7 kl = 

(iii) The trace in the graphical calculus is given by appropriately joining wires to close 
loops. 



Together the cups and caps give rise to a correspondence between different types of maps 
and states. We call the duality induced by bending and exchanging wires Penrose duality. 
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2. PENROSE WIRE BENDING DUALITY 



Now we will consider the set of operations formed from bending tensor wires forwards or 
backwards using cups and caps, as well as exchanging wires using SWAP. We can conceptu- 
alize this set of transforms acting on a tensor as amounting essentially to matrix reshapes. 
From the snake equation, action with a cup or a cap is invertible and SWAP is self in- 
verse. This implies that all possible configurations of a tensor's wires obtained using these 
operations can be considered equivalent. We will start with an example. 

Example 1. Given a tensor T 1 - with fixed labels i,j one uses cups and caps to rearrange 
index elevations, arriving at 

r,. v". t„. •/;.•' (5) 

Using the SWAP operation one reorders the horizontal position of i and j. Then applying 
cups and caps yields 



rp] rp'jl rp rp 



(6) 



for a total of eight possible reshapes. 

For an n index tensor, each index can be either up or down, yielding 2 n possibilities. The 
symmetry group formed by SWAP is of order n\ and acts to arrange the horizontal position 
of the n legs of a tensor, yielding (provided we distinguish forms of the type T { 1 and T\ in 
Example [T]) n! • 2 n different ways to reorder the n indices of a tensor. 

Remark 2 (Ordering operators by numbers of inputs and outputs). In the previous example, 
we considered T i 1 (b) and T\ (a) as distinct. This is illustrated in (a) and (b) as follows. 



(a) 



(b) 



This provides an example of an awkward property of the inherently one-dimensional Dirac 
notation. Both (a) and (b) represent the same map, but when we write them in a basis, 
consistency dictates that one will expand in the basis (i\ <S> \j). 

With this equivalence in mind, we note that the tensor T ^ from Example factually has six 
unique reshapes, as two pairs of reshapes are diagrammatically equivalent. In other words, 



in (b) below we have that T l - 
(a) 



(d) 




rpi rp 

3 ~ 3 
(b) - 



and for (e) we have the equality T- 
— T„ l (c) 



T{ = T, 



(e) 



T 



(11 




rp% 3 



Together, we call these reshapes the natural tensor symmetry class. More generally one 
finds the number of diagrammatically unique reshapes of a tensor by (i) counting the number 
of possible ways it can have its wires bent, either forward or backwards using the cups and 
caps, and (ii) the number of ways a tensor can have its arms and/or legs reordered. We 
arrive at the following result: 

Theorem 3 (Natural tensor symmetry class) . The arms and legs of a tensor r*£.'.'.« with n 
input and output legs in total can be rearranged in (n + 1)! different ways. 
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3. DIAGRAMMATIC SVD 



In this section, we introduce the diagrammatic representation of the singular value de- 
composition (SVD). Later it will be used to simplify invariants obtained through network 
contraction, and iterated to obtain a matrix product state (MPS) description for a pure 
state. 

The SVD factors tensors into well defined building blocks with simplistic interaction 
properties: (i) a valence-one tensor storing singular values, (ii) a valence-three-COPY tensor 
used to create a diagonal map, and (iii) a pair of valence-two unitary gates. COPY-tensors 
have been studied in the setting of the Penrose tensor calculus, in work dating back at least 
to Lafont 



10, 11 



see also [12, LL5L 16 



Definition 4 (COPY tensor). The m-to-n COPY tensor is defined in the computational 
basis as 



d-i 



COPY m ^ n := £ | k_j_k)(k_j_k 



(7) 



k=0 



It is named accordingly because connecting a basis state \k) to any of its input or output 
wires collapses the sum and breaks the tensor up into unconnected copies of \k) and (k\. As 
with classical circuits, in the diagrammatic tensor notation a COPY tensor is represented 
by a simple black dot • with m input and n output legs. For a brief enumeration of its 
algebraic properties, see [151 . 



Theorem 5 (Diagrammatic SVD). Any valence-two tensor f : A —t B can be factored 
into a non-negative, unique valence-one tensor E ; an valence-three COPY tensor, unitary 
valence-two tensors U and V , and a diagonal valence-two dimension changer tensor Q when 
necessary, as shown in Figure® 




FIG. 2. Diagrammatic singular value decomposition for / : A — >■ B. The dimension changing 
tensor Q is a rectangular tensor which has all l's down the diagonal and zero entries otherwise. 



Proof. The SVD of / is 

/ = uzv, 

where U : B —> B and V : A — >■ A are unitary operators and E : A — > B is diagonal in the 
computational basis, with the (necessarily non-negative) singular values cr^ of / along the 
diagonal. E can be written as 

d-1 d-1 
S = ^2 a j\j) B (j\A = ^2\ i )B( i \A^2\j)AUj\A^2 (T k\k) A (<T k > 0), 




Qab COPY 2 _>i 
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where d = min(dim A, dim 5). We have expressed £ as a contraction of an valence-one 
tensor a with a COPY tensor. The non-square tensor Qab is only necessary if A and B 
have different dimensions. □ 

Corollary 6 (Diagrammatic Schmidt decomposition). Given a bipartite state we use 
the snake equation to convert it into a linear map f : A — >■ B ( inside of the dashed region 
in Figure [fl Now we apply the SVD as in Theorem^ Diagram reorganization leads to 







Q 



B 



Uhg- 



FIG. 3. Diagrammatic Schmidt decomposition for \tp) £ A® B. 



the diagrammatic Schmidt decomposition of 
to the Schmidt coefficients. 



The singular values in a now correspond 



The network topology of the diagrammatic Schmidt decomposition can be used to study 
the entanglement properties of the bipartite state 

Example 7 (Entanglement topology). The topology of a bipartite state depend on the 
singular values in the triangular tensor \a) = (Tq\0) + &\ | 1) + . . . + cr^^d — 1). The most 
significant topology change occurs when the input state to the black COPY tensor is single 
basis state \a) = \k) — this causes the diagram to break into two halves, each storing a 
copy of \k) as shown in Figure 11(c). When the input state is a unit for the COPY-tensor, 
the tensor structure is converted to a smooth wire (this is the maximally entangled case), 
as illustrated in Figure 0|b). 



(a) J— I (b) J-"L 

«R cR 



L (c) 



40 



FIG. 4. Tensor network entanglement topology, (a) the general form of a bipartite state \ip) = 
0i\(pi)\(f)i} . The state takes the form (b) iff the Schmidt coefficients all take the same value. In 
this case, the state can be transformed using local unitary operators to the generalized Bell state 
in dimension d. The state takes the form in (c) iff only the first singular value is one. In such a 
case the state is separable, as evidenced by the diagram breaking into two pieces by applying the 
rewrite rules of the COPY-tensor. 



2 This has also been understood as a diagrammatic form of map-state duality underlying bipartite entan- 
glement evolution [ijj]. 
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4. DIAGRAMMATIC MATRIX PRODUCT STATES 



As a key application of the diagrammatic SVD, we will consider matrix product states (MPS), 
an iterative method to factor quantum states into a linear chain of tensors (see 0, 0, Hl[ 
for a recent review and [l5[ for work considering the category theory behind MPS). The 
reason this factorization is called a ID method is because it is known to describe a class of 
ID systems efficiently, and because the factorization results in a ID chain (for a discussion 
of other factorizations and the connection to geometry see jll[). Without loss of generality, 
we will apply the MPS method to a four-party state, and explain the procedure in terms 
of three distinct steps. Consider a quantum state, expressed as a triangle in the Penrose 
graphical notation with a label 1 inside and open legs labeled k,m. 




i j k m 



(Step I). We will create a bipartition comprising a first collection, containing only leg 
% and a second, containing legs j, k, m. We will then apply the diagrammatic SVD across 
this partition. The partition is illustrated with the dashed cut below in (a). Figure (b) 
results from applying the diagrammatic SVD across this partition, factoring the original state 
labeled 1 in (a) into a valence- two unitary box labeled 2, a valence-one triangle containing 
the singular values labeled 3, and a valence-four triangle labeled 4, all contracted with a 
COPY-tensor, as illustrated. A new internal label (d) for the wire connecting the COPY- 
tensor to the valence-four triangle (4) was introduced for clarity, (see also Figure S£ a) and 

hd>))- 




Remark 8 (Isometric internal tensors). The valence- four triangle tensor in (b) above arises 
from contracting a unitary map with a dimension changing tensor Q (see (a) below). The 
input leg shown is labeled d. The other legs are contracted with a fixed basis state |0), from 
the SVD in (a) above. We then depict this as the triangle 4 as in (b) above. From the 
unitarity property, the isometry property follows, as illustrated graphically in (c) and (d) 
below. 
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(a) 





j k m 



j k m 



(d) |A 




B — 



(Step II). To illustrate the next step in the factorization, we will remove the tensor labeled 
4 by breaking the wire connecting it to the COPY-tensor (a). We will then partition this 
separate tensor into two halves, one containing wires d, j the other half wires k,m. This 
partition is illustrated by placing a dashed line (labeled cut) in (a). We arrive at the the 
structure in (b), which we have explained in the first step, (see also Figure [3Tb) andHIc)). 



(a) 




(b) d 



T 



W/i 
1 



k m 



Remark 9 (An elementary property of tensor network manipulation). It is a fundamental 
property of tensor network theory that one can remove a portion of a network, alter this 
removed portion of the network without changing its function, and replace it back into the 
original network, leaving the function of the original network intact. 

(Step III). In the third and final step of the MPS factorization applied to this four-party 
example, following remark [9] we first place the tensor we have factored in the second step, 
back into the original network from the first step, see (a) below. We then repeat the second 
step, applied to the triangular isometry tensor, labeled internally with a 7. This results in 
the factorization appearing in (b). (see also Figures HJc) and@Jd)). 

(a) cut (b) 



j k ; m 



w TvA " w Tw 




Remark 10 (Step n). The iterative method continues in the same fashion as the first three 
steps, resulting in a factorization of an n-party state. A summary of the MPS factorization 
applied to a four-party state is shown in Figure HI 

(Summary). We will now consider Figure HJ which summarizes the factorization scheme. 
In the steps we have outlined, we have factored Figure HJa) into the MPS in Figure 11(d), in 
terms of the components listed below. 



(i) States (labeled 3, 6 and 9; denoted 3 , c, 
06 = (A 2 , A 3 , A 4 , A 5 ) T and </> 9 = (A 6 ,A 7 ) T . 



i and (f)g, respectively): </> 3 = (Ao,Ai) T , 
The Aj's are the singular values across 
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cut 



(b) 



cut 

d i (<0 



cut (d) 




j k m 



T7 



w \w 



k m 







i J 


k 



FIG. 5. MPS factorization steps (diagrammatic description of steps I, II and III). The quantum 
state (a) is iteratively factored into the ID matrix product state (d). This procedure readily extends 
to n-body states. 



each partition. The number of non-zero singular values (x) is given by the minimum 
dimension of the two halves from the cut. For the case of qubits, the first outside 
partition has at most two non-zero entries, and the next inside partition has at most 
4. One might also consider the singular values as the square roots of the eigenvalues 
of either member of the pair of reduced density matrices found from tracing out either 
half of a partition. 

(ii) Unitary gates (labeled 2 and 10; denoted U2 and Uiq, respectively). 

(iii) Isometries (labeled 5 and 8; denoted I5 and /§ respectively). The isometry condition 

describes the tensor relation Ij I J r q = 5 d r . It is a consequence of the fact that tensors 
5 and 8 arise from unitary gates, as explained in Remark [HJ The isometry condition 
plays a more relevant role in structures other than ID tensor chains. 

We note that by appropriately combining neighboring tensors as in Figure 0Ja), one 
recovers the familiar matrix product representation of quantum states BJb). Matrix product 
states are written in equational form as 

* = E A t ] Af A k ] A [ £\ijkm) (8) 

i,j,k,m 

Here A^ becomes a new tensor formed from the contraction of tensors labeled 2, 3, and 
is a contraction of tensors labeled 5 and 6, etc. The exact grouping choice has some 
ambiguity. 

A utility of our approach summarized in Figure 0Ja) is that the COPY-tensor is well 
defined in terms of purely graphical rewrite identities (as seen in Definition [4]). These 
graphical relations allow one to gain insights (into e.g. polynomial invariants as will be 
seen). The factorization we present however, allows one to preform many diagrammatic 
manipulations with ease, and exposes more structure inherent in a MPS. 

Remark 11 (Data compression). The compact representation of a MPS is recovered by 
picking a cutoff value for the singular values across each partition, or a maximum number of 
allowed singular values. This allows one to compress data by truncating the Hilbert space 
and is at the heart of MPS computer algorithms in current use. 

The singular values found from the MPS factorization can be used to form a complete 
polynomial basis to express invariant quantities related to an MPS. 
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(a) 



AS- 

h5 



mrwm 



■a 



(b) 



AM 



FIG. 6. Conversion from our notation (a), to conventional MPS notation (b). The factorization 
methods we have presented here and elsewhere la, [l6| allow one to "zoom in" and expose in- 
ternal degree of freedom (a) or "zoom out" and expose high-level structure (b). The equational 
representation of the MPS in (b) is given in ©. 



PENROSE NOTATION MEETS ENTANGLEMENT INVARIANTS 



Here we will consider the variant of the graphical tensor calculus of Penrose [6|] we have 
tailored to represent and contract polynomial invariants. We must first recall the notions 
surrounding polynomial invariants. 



Polynomial invariants 

Assume we are given a group G, a vector space V, and a group representation D : G — > 
Given a set Q, a function / : V — > Q is an invariant function or simply an 
invariant under D iff it is constant on the orbits of 

In the context of quantum mechanics, the vector space V is typically either the state 
space H, or End("H), the space of linear operators H —>7i. Any representation D : G —> Aut(H) 
on % induces a representation R : G — > Aut(End("H)) on End("H): 

R(g)(p) = D(g)pD~\g). (9) 

An important class of invariants are the polynomial invariants / : V — > C, which are 
polynomial functions of the coefficients of p or in the standard basis. The study of such 
polynomials is known as invariant theory [22[ . David Hilbert made notable progress on this 
topic, which he perused throughout his life. There has been past work considering these 
invariants in the context of quantum information science. Some that was influential to us 
includes 

Remark 12 (Basis independence). To form a polynomial out of the coefficients of a state, 
one first chooses a basis to express the state in. The value of the polynomial generally 
depends on the basis chosen. However, a polynomial that is invariant under any group that 
contains the local unitary group as a subgroup is inherently basis independent as long as 
our basis is a tensor product of orthonormal local bases. 



Invariance under the local unitary group 

In this section we will study local unitary (LU) invariants. 



3 Aut(V) denotes the group of automorphisms of V, i.e. the invertiblc linear maps from V to itself. 

4 Or equivalently iff / itself is a fixed point under the induced representation D' : G — > Aut(F(V, Q)), 
D>(g)(f) = foD(g-i). 
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Definition 13 (LU equivalence of states). Two quantum states (pure or mixed) in the 
Hilbert space % — %i ® %2 ® ■ ■ ■ ® 7~L n are LU equivalent iff they are related by a local 
unitary transformation, that is, a member of the natural representation of the group 



G 



LU 



U(l) x SU(dx) x SU(d 2 ) x ... x SU(d n ), 



(10) 



where di = dim'Hj is the dimension of the ith subsystem. LU equivalence yields a par- 
titioning of the state space into LU orbits. Entanglement measures are by definition LU 
invariants, i.e., constant on the aforementioned equivalence classes. 

We now present a diagrammatic method for systematically generating polynomial LU 
invariants for state vectors and operators by casting the method of Grassl et al. [l|, |23j into a 



form based on the Penrose tensor calculus. The method generates homogeneous polynomials 
in the state coefficients that are necessarily invariants of the local unitary groupH A utility 
of generating this set stems from the fact that the tensor networks considered can be used 
to calculate quantities that are invariant under the action of the local unitary group. 
Given a density operator p : 7i — > Ti, consider the network 



equivalent to the expression Tr(Fp) = F l jP > i - By choosing a suitable F, we can represent 
any first-degree homogeneous polynomial in the coefficients of p in this way. Likewise, the 
tensor network 



P - 



translates to Tr(Fp® 2 ) = F d j k p > i p k l , giving all the second-degree homogeneous polynomials. 
The procedure carries on in this fashion. A diagram with k copies of p gives us all the 
homogeneous polynomials of degree k: 



Tr(Fp 5 



pq-t 



p P ip q j ■■■P t 



(11) 



Having thus generated a complete basis for the polynomials in the coefficients of p, we next 
wish to find out which of these homogeneous polynomials are invariant under the natural 
representation of G\jj: 



Tr(F(UpU- l ) m ) 
This is fulfilled iff 



Tr((U- 1 )^ >k FU m p 



Tr(Fp 



VUeG LV , V P . 



;i2) 



[F,£/® fc ] = V[/GC7 LU . (13) 

We are then faced with finding matrices F that commute with U® k for each U G G*lu- The 
solution is roughly stated in the following theorem. 



5 Although we can generate a complete set of invariants in this fashion, except in rare cases, finding the 
minimal complete set of polynomial invariants is computationally difficult. This alternative line of research 
has been a key focus in the connection of invariant theory with quantum entanglement. 
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Theorem 14 (Brauer |24j, Procesi [23]). The algebra of matrices that commute with every 
jj®k j Qr jj ^ Glu is generated by the unitary representation 

T: S k x ... x S k -> Aut^ fe 

S v ' 

n copies 

of the n-fold direct product of the permutation group S k which, independently for each of 
the n subsystems, permutes the relative ordering of the k copies of that subsystem's state 
space within the total space 1-L® k . 

Hence, it is enough to consider matrices F which correspond to these permutation maps. 
The permutation group has a well known and evident diagrammatic form. Below, we show 
the elements of the permutation group (a) Si, (b) S2, and (c) S3. 



(a) 



(b) 



X 



(c) 




X 




We then carry on to evaluate all the expressions of the form 

h;a u a 2 ,...,a n (p) ■= Tr(T(ai, cr 2 , . . . , a n )p® k ), where cr, <G S k (14) 

to generate the homogeneous invariant polynomials of degree It is easy to see that the 
diagrams we obtain are indeed invariant under Glu ; as shown in Figure 



(a) 



(b) 



G 



G 



u 


p 


w 


u 


p 


w 




u-u 



u-u 



J 




(c) 



= GrrTrrU 



FIG. 7. Proof of the invariance of i2;(i2)- Having acted on p with some unitary operation U, we 
slide U and W around the bends, taking the transpose and resulting in (b). The unitaries cancel 
and the diagram reduces to (c), showing that it indeed describes an invariant. A little bit of further 
manipulation shows that ^2;(i2) evaluates to Tr(p 2 ). 

Not all of these invariants are independent, or even distinct. We can eliminate some of 
the redundancy using the following theorem, proven in [l[: 

Theorem 15 (Invariant distinctness Since all the copies of p in (fT4|) are identical, 
we may permute their relative order without changing the invariant. This is equivalent to 
conjugating each subsystem permutation G{ with the same element r e Sk-' 



k'TcriT 1 ,...,T(T 7l r 



V<Tj,r G Sk. 



(15) 



This theorem enables us to arrange each invariant diagram to the following canonical 
form which makes it easy to tell if two diagrams are topologically distinct. 



We use the cycle notation to denote specific elements a of the permutation groups. 
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5 PENROSE NOTATION MEETS ENTANGLEMENT INVARIANTS 



(i) The k copies of the system are arranged such that the permutation on the first sub- 
system is grouped by cycles, ordered by non-increasing cycle length. 

(ii) The process is repeated on the second, then third etc. subsystem within the remaining 
permutational freedom, i.e. cyclic permutation within the cycles and permuting cycles 
of identical length. 

If a particular diagram is not connected, the corresponding invariant is the product of the 
invariants corresponding to the disjoint subdiagrams. 

Remark 16. Note that using the procedure here, two algebraically independent invariants 
necessarily have topologically distinct diagrams but the converse does not necessarily hold. 

Theorem 17 (Real-valuedness of the invariants). If all the permutations Oi are self-inverse, 
or can all be inverted by conjugating them with the same element r as shown in Theorem \Tb\ 
the invariant Ik-ai,...^ i> s necessarily real for all states. 



Proof. 



It,^..,*M =Tr((^*)trt( ffll ... |fffl )) ^Tr^t^,...,^)^ 
= Tr (T(cr 1 _1 , . . . , V*) = W „.»• 



(16) 



□ 



Theorem 18 (States with a single subsystem). All the independent invariants of a d- 
dimensional state p with a single subsystem are given by 



4:=Tr(p fc ), where k e {1, 2, . . . , d}. 



(17) 



Proof. The only degree one invariant, I\ := Tr(p), is presented in Figure El a. The two 
possible degree two invariants are shown in FigureElb. The first one is simply If. The second 
one, I2 '■= Tr(p 2 ), however is independent. Likewise, the only independent invariant of 
degree three, ^3 := Tr(p 3 ), is given in Figure EJc. In general, at each degree k we obtain a 
single new independent invariant Ik '■= Tr(p k ), by using a complete permutation connecting 
all the k copies of the state. 



(a) 



(b) 



CnrP G 



(_z_) G 



(c) 





FIG. 8. LU invariants for a single subsystem, (a) I± := I\. e = Tr(p). (b) 1^ 



II h := I 



2;(12) 



Tr(^). (c) J 3 :=J : 



3;(123) 



Tr(^). 



The Cayley-Hamilton theorem now tells us that the basis is finitely generated, as every p 
satisfies its own characteristic polynomial, giving an <ith degree polynomial equation in p 
which enables us to express any I m with m > d in terms of the lower-degree invariants [23J . 

□ 
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5 PENROSE NOTATION MEETS ENTANGLEMENT INVARIANTS 



Example 19 (Invariants for a single qubit). The only independent (fundamental) invariants 
of a single qubit state are ii and I2, defined in the previous theorem. I\ = Tr(p) is the norm 
of the state. I 2 = Tr(p 2 ) turns out to be precisely the purity of the state. In terms of the 
eigenvalues (Ao, Ai) of p we have I\ — Ao + Ai = 1 (for normalized states) and I 2 = A 2 , + A 2 . 
From the Cayley-Hamilton theorem we have that there is a second degree monic polynomial 
in p that vanishes identically. In other words, constants a, b exist such that 

p 2 + ap + 61 = 0. (18) 

Multiplying both sides by p m and taking the trace, we obtain the recurrence relation 

Im+2 + al m+1 + bl m = 0, (19) 



and thus find that the traces of higher powers of p can be expressed in terms of I\ and 
I 2 . These invariants are indeed algebraically independent and complete, meaning any other 
polynomial invariant can be expressed in {R, +, •, I\, I2}. For instance, 



det(p) = \ (Tr(p) 2 - Tr(p 2 )) = ~ [l\ - I 2 ) = AqAl 



Likewise, J 3 = Aq + A^ can be written as 



h(I 2 -det(p)) 



(20) 



(21) 



Remark 20 (Bipartite states). For bipartite states we obtain a much more complicated set 
of invariants. Figure M presents all the topologically distinct invariants up to k = 3. 



(c) C 



(b) a 



b 



D 







G 



^ G 



3 



G 
d 



FIG. 9. LU invariants of a bipartite system up to k = 3. To avoid listing essentially similar diagrams 
we only show here the distinct invariants modulo swapping the order of the two subsystems, (a) The 
only first degree invariant is Ii- e ,e = Tr(p), same as with a single subsystem, (b) In the second 
degree we obtain a new invariant, i2;(i2),e- ( c ) There are several new third degree invariants, 
including the topologically distinct / 3; (i23),(i23) and 4 ; (i23),(32i) • 



Pure states 



If the state p is pure, the diagrammatic structure of the LU invariants simplifies consid- 
erably, and many of the diagrams break up into unconnected subdiagrams. Furthermore, in 
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5 PENROSE NOTATION MEETS ENTANGLEMENT INVARIANTS 



the case of bipartite pure states, we may apply the Schmidt decomposition and introduce 
graphical rewrite rules to show that these invariants reduce to polynomials of the Schmidt 
coefficients. 

Theorem 21 (Bipartite pure states). Applying the diagrammatic Schmidt decomposition 
presented in Figure d to the bipartite invariant diagrams in Figure 0, we can see that the 
unitaries U , V and the dimension changers Q always cancel, and the invariant diagrams 
break up into mutually disjoint loops corresponding to sums of even powers of the Schmidt 
coefficients {o~i}flo- Hence, the only distinct invariants we obtain are of the form 

Jk '■= h; (12-fc),e = ^2 a i k ' 
i 

for /r6 {1,2,...}. In Figure [TOl we present this process for the invariant 




= <^S2) = °C 

FIG. 10. Invariant I3 ; (i23),(i2) evaluated for a pure bipartite state p = \ip){ip\ using the diagram- 
matic Schmidt decomposition. The unitaries U and V (and possible dimension changers Q) cancel, 
and one is left with two disjoint loops, on which the blank circles denote diagonal tensors with the 
Schmidt coefficients {ai}f~Q on the diagonal. 

Since the d Schmidt coefficients themselves (by construction) form a complete set of 
bipartite LU invariants, we should be able to express them as functions of {Jj}f =1 . This is 
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5 PENROSE NOTATION MEETS ENTANGLEMENT INVARIANTS 



accomplished in principle by solving the following system of polynomial equations: 

i 



i 

Example 22 (Pure state of two qubits) . Given a pure two-qubit state 



(23) 



we wish to compute its LU invariants. Since it has two Schmidt coefficients, we expect to 
find two independent invariants. The first one, J 1; corresponds to the squared norm of the 
state, and is thus trivially invariant under unitary transformations of \ip): 



The second independent invariant, J 2 , doesn't have as simple an interpretation: 



J 2 = 




The diagrams correspond to the equations 

J x = Y^aFaiij = al + of, 



with the solution 



ijkl 



<i=2 ( Ji±\/2J 2 -J 2 1 ), 



(24) 



which yields J 2 < Jf < 2J 2 - For a normalized state {J\ = 1), it can be shown that the 
invariant J 2 can be expressed as 



h — 1 - 2|a ottn - a io;io| 



(25) 



where aoo«n — ttoi«io is simply the determinant of the coefficient matrix a, which is non-zero 
iff represents an entangled state. 
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5 PENROSE NOTATION MEETS ENTANGLEMENT INVARIANTS 



Ji and J2 are the only algebraically independent LU invariants of a pure two-qubit system. 
Any polynomial function of such invariants is also a polynomial invariant. In this fashion, 
it is a remarkable feature that functions of J\ and J2 are all that is needed to express any 
local unitary invariant of two-qubit pure states. This elementary result follows from a much 
more powerful and general result in classical invariant theory, a proof by Hilbert that the 
ring of polynomial invariants is finitely generated [22|. This corresponds to freely generated 
linear sums and products of J±, J2, e.g. the ring {Ji, J 2 , (K, +, •)}■ Any minimal complete 
set of invariants that can freely generate the full ring are called fundamental invariants. 



INVARIANTS, ENTROPIES AND ENTANGLEMENT 



Here we focus on expressing Renyi entropies in terms of the invariants we have found 
tensor contractions for in the previous sections. The Renyi entropy has many uses in con- 
densed matter physics (see for example [HI, [26| ) and recently has been given an interesting 
physical interpretation j27| . We will recall the definition as 



Definition 23 (Renyi entropy [28] ). The Renyi entropy of order a is defined as 

1 



1 



Q 



lnTr(p Q ). 



In the limit a — > 1 we obtain 



lim S a = — Tr(pln p), 



a— >1 



(26) 



(27) 



which is the von Neumann entropy. 



Here we note that terms such as Tr(p°) are in correspondence with the tensor contractions 
evaluating to invariants which we have already found. To explain how we can contract tensor 
networks to evaluate Renyi entropies for counting a > 1, we will close the paper with a 
specific example, though the procedure we describe is general. We will focus explicitly on 
the invariants of a bipartition of a 5-party qubit state \ip) G C 2 ® C 2 ® C 2 ® C 2 ® C 2 . We 
first recall that we can factor any state into a MPS; in the case of our example, this yields 
the following graphical depiction. 



(a) 



(b) 



1 1 1 1 1 



Here (a) is the original state, (b) is our factorization in terms of COPY-tensors, and (c) 
recovers the familiar MPS representation, as explained in Section |H The method to evaluate 
Renyi entropies by tensor contraction works generally by first grouping the legs of any tensor 
network state into a bipartition, and to consider the correlations between the two halves. In 
the present example, we group the two top legs (^4) and the other three legs (B) and then 
apply the graphical SVD: 






Q- 



- B 
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After tracing out system B, we are concerned with a four dimensional space C 2 ® C 2 . From 
the Cayley-Hamilton theorem 

Tr(p 4 ) + aTr(p 3 ) + 6Tr(p 2 ) + c = (28) 

for some values of a, b, c. We hence conclude that all information we can expect to find 
can be obtained by evaluating tensor contractions for Tr(p n ) for n = 2,3,4 as given in the 
section on invariants. As an example, to evaluate Tr(p 3 ): 



We can evaluate this contraction using a tensor network numerical algorithms package. By 
writing it in terms of the graphical SVD and applying the rewrite rules developed in the 
present work, we arrive at 




This illustrates that the network reduces to an expression in terms of the singular values 
of the pure state. If the singular values of the original state are given as (a±, cr 2 , cr 3 , a 4 ) =: 
(y/P~v Vp~2> V^3' y/P/^i then the expression evaluates to 

/ 3 ;(i23), e = a\ + 4 + 4 + °t=p\+Pl+Pl+Pl (29) 

where the eigenvalues Pi on the right are the probabilities of measuring the reduced system 
in the ith eigenstate. We then express the Renyi entropy with a = 3 as £3 = — |ln(/3 ; ( 123 ) ie ). 
Other quantities can be similarly calculated, resulting in the following identical relation for 
the Renyi entropies 

exp(-35 4 ) + aexp(-2S 3 ) + bexp(-S 2 ) + c = (30) 
which is an alternative expression for (T28]) in terms of the S a from Definition [23j 

6. CONCLUSION 

We have developed tensor contractions to express a complete polynomial basis for the local 
unitary invariants of any finite-dimensional quantum system. Using the diagrammatic SVD, 
we proved that for pure bipartite systems these contractions can be expressed in terms of 
manifestly invariant singular values. These methods seem to provide new conceptual insight 
to understand quantities such as entropies which are expressed as a holomorphic function of 
the invariants. By connecting invariant theory with tensor network states, we hope that this 
work leads to a better understanding of how entropies and entanglement measures can be 
calculated for specific models, and to ultimately lead to a better understanding of invariants 
in various higher-dimensional geometries as well. 
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